Half-integer contributions to the quantum Hall conductivity from single Dirac cones 
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While the quantum Hall effect in graphene has been regarded as a realization of the anomaly 
associated with the massless Dirac particle carrying half the usual topological integer, this is hidden 
due to the doubling of the Dirac cones. In order to confirm the half-integer contribution from each 
Dirac cone, here we theoretically consider a lattice model in which the relative energy between the 
two Dirac points is systematically shifted. With an explicit calculation of the topological (Chern) 
number, we have demonstrated that each Dirac cone does indeed contribute to the Hall conductivity 
as the half odd integer series (• •• , —3/2, —1/2, 1/2, 3/2, • ■ ■ ) when the Fermi energy traverses the 
(shifted sets of) Landau levels. The picture is also endorsed, via the bulk-edge correspondence, from 
the edge mode spectrum for the present model. 
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I. INTRODUCTION 

In the seminal discovery of the quantum Hall effect 
(QHE) in graphene 1 - 3 , a most striking point is that 
the graphene QHE is regarded as a realization of the 
anomah^ in the massless Dirac particle, where each 
Dirac cone carries 1/2 of the usual QHE topological 
number 5,6 . Namely, in place of the usual Hall conduc- 
tivity, a xy — 0, 1,2, ••• (in units of —e 2 /h, with the 
spin degrees of freedom dropped here), we have a xy — 
2n + l(ra = 0, ±1, • • • ). The standard understanding is 
that the honeycomb lattice has a couple of Dirac cones at 
K and K' points in the Brillouin zone, so that we just di- 
vide the above formula by 2 to have a xy = n+1/2 for the 
contribution from each valley. The half integers are thus 
hidden in the total Hall conductivity. Since the Nielsen- 
Ninomiya theorem! dictates that Dirac cones should al- 
ways appear in pairs as far as the chiral symmetry is 
present, we can pose an important question: is it really 
impossible to resolve the half-integer components? 

In terms of field theory (as opposed to lattice mod- 
els), the situation is simple: the Hall conductivity for 
a massive (m 7^ 0) Dirac particle, in zero magnetic 
field, is given by <r xy — (l/2)sgn(m) when the Fermi 
energy Ep is in the mass gap 4 . This can be readily 
shown by calculating Berry's connection. If we con- 
sider a field theoretical model concerning a single Dirac 
fermion at kg with a gap \m\, the Hamiltonian is given 
by h = 5k x a x + 6k y a y + ma z = R ■ cr, where R = 
(i?i , i?2 , R3 ) = (Sk x ,5k y ,m), 5k = k — k^. If we denote 
the eigenstate having an energy —R= — ^\5k\ 2 + m 2 
as tp, and its Berry connection as A = -ip^dip, the total 
Berry curvature over the (k x ,k y ) plane (represented as 
n m ee {R\R 3 = m}) is C = ±i / nm dA = f Um B ■ dS = 
sgn (m)/2, where B = Rj (47ri? 3 ), which is the flux from 
a magnetic monopole at the origin with div_B = S 3 (R) 
and dS = dk x dk y dk x R X dk y R = dk x dk y (0,0,l). We 
can then see that we have 1/2 since just the half the 
total flux of the magnetic monopole passes through the 



plane Tl m - This implies that the topological change in the 
quantum ground state (filled Dirac sea) is characterized 
by Aa xy — ±l 9 ~ n . Generically speaking, a topological 
quantum phase transition is naturally accompanied by a 
sign change in the mass in the effective, low-energy Dirac 
fcrmions. Quantum Hall plateau transition is a typical 
example, and the topological insulators such as the quan- 
tum spin Hall system are also described along this line, 
where the spin-orbit interaction induces the sign change 
in the mass 1 ^— . 

While the quantization into the half odd integers 
is conceptually interesting, if we go over to lattice 
models, however, we have a periodicity in the Bril- 
louin zone, which implies that the topological numbers 
should always be integers as dictated by the celebrated 
but inescapable Thouless-Kohmoto-Nightingale-den Nijs 
(TKNN) formula 5 -. The integer total Hall conductivity 
for graphene with a pair of Dirac cones is an example of 
this. If we turn to a wider class of lattice models, we can 
go around the Nielsen-Ninomiya theorem. For instance, 
we can have a lattice model that has an odd number of 
massless Dirac cones, but even in that case we still end up 
with integer Hall conductivities, which may be regarded 
as due to hidden massive Dirac fermions required to exist 
for the topological consistency in the lattice system to be 
guaranteed. The massive Dirac fermions (sometimes re- 
garded as "spectators"—) are hidden in the high-energy 
region, which do not appear in the low-energy physics 
except that they make sure that the the Hall conductiv- 
ity are topologically protected to be integers. In another 
manipulation of Dirac cones^, a quantum Hall effect in 
zero total magnetic field has been considered, where the 
Hall conductivity is shown to take the value of e 2 /h even 
in zero field in a model containing complex hoppings in 
the situation where the Dirac cones are made massive. It 
is thus rather difficult to confirm half-odd-integers in the 
Hall conductivity for a single Dirac cone. 

In this Rapid Communication, we want to shed light 
to this problem by posing the following question: if we 
can manipulate the energies of the multiple Dirac points, 
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FIG. 1: (Color online) (a)The honeycomb lattice model con- 
sidered here with nearest-neighbor hopping to (solid lines) 
and second-neighbor hopping which is iti along each dashed 
arrow, —iti in the opposite direction. A unit cell and A, B 
sublattices are indicated. (b)The dispersion relation in the 
present model, seen horizontally, with Dirac cones shifted in 
energy by A = v^ti . 



can the half integer series confirmed through a system- 
atic behavior of the total Hall conductivity? Namely, 
we shall construct a lattice model, where the two Dirac 
points arc shifted with the massless cones preserved for 
both of the Dirac points. We shall identify, by di- 
rectly computing the topological Chern number for sys- 
tematically varied relative position of the Dirac cones, 
that each Dirac cone indeed has a half-odd-integer series 
(••• ,-3/2, -1/2, 1/2, 3/2, •••) when E F traverses (now 
shifted sets of) Landau levels belonging to the two Dirac 
cones. The picture is further endorsed, via the bulk-edge 
correspondence, from the behavior of the edge modes for 
the shifted Dirac cones for finite systems. 



II. SHIFTED DIRAC CONES 

A strategy to conceive a model in which two Dirac 
cones are preserved but mutually shifted in energy is the 
following. A simplest solution is to add a term that is 
proportional to <jq (unit matrix) in the space of Pauli 
matrices with a fe-dependent coefficient. Thus we can 
introduce a Hamiltonian, 

H = J2^2c{ a [hf + 2t 1 (smk 1 )a a } a ^c kf 3, (1) 

k o,j8 

h k r = t [(1 + cos ki + cos fc 2 )i7i + (sin k\ + sin £2)02] , 



where c k creates an electron in k- space, and <7j's are 
Pauli matrices with a, /3 denoting their components. In 
this Hamiltonian, we have added, on top of the nearest- 
neighbor hopping tg, an extra exo term with hopping t\ 
as a coefficient. This only affects the band dispersion 
as an additive, fe-dependent term, which implies that we 
lift the degeneracy between K and K' points if the fe- 



dependent term has different values at K and K'. A sim- 
plest choice is oc sinfci. If we go back to the real space, 
the tight-binding model is as depicted in Fig |l(a)| which 
has extra second-neighbor hoppings. The added hopping 
has to be only between A- A and B-B for the Dirac cone 
to be preserved, and they have to be complex for the 
degeneracy between K and K' (mutually time-reversal 
pairs) to be lifted. So the model with a complex hopping 
is rather fictitious, but we do accomplish shifted Dirac 
cones as depicted in Fig 1(b) In this model the chiral 
symmetry [i.e., the presence of an operator (a z in the 
case of graphene) that anticommutes with the Hamilto- 
nian] is broken, since appearance of o§ in % invalidates 
{W,a z } = 0. Nevertheless, the addition of ctq preserves 
the shape of Dirac cones, along with the species dou- 
bling. If we expand Hamiltonian ([1]) around fco (K or K' 
points), we have 



hi 



(2) 



where \ = +1(— 1) correspond to K (K'), A = \/Zt\ is 
(half) the shift, vp = V3atQ/2h the Fermi velocity, and 
5k = k — fco. So the effective theory around these points 
are a pair of Dirac equations at K and K' with a shift 
±A in energy as desired. 



III. CHERN NUMBERS 

The Landau levels in a magnetic field B z expected from 
the effective Dirac Eq. ([2]) are = \IS. + ?iaj c sgn(n)- v /n 
[ui c = 3^H ^/2^, <j) = {V3a 2 B z /2)/{hc/e),n = 
0,±1,±2,---]. However, we should of course go back to 
the original lattice model for obtaining the Landau levels 
and, especially, the Hall conductivity. The calculation 
of the Hall conductivity, which is a topological Chern 
number, requires a bit of care, since we have to question 
the behavior of the topological number around E = 
which is just where the hole and electron branches meet 
in the Dirac cone dispersions. This implies we have to 
sum over all the contributions from the "Dirac sea" to 
get the total Hall conductivity. This ordinarily poses a 
numerically difficult problem, but we can overcome the 
problem with a method that employs a non-commutative 
Berry's connection^ and its integration (Chern number) 
over the Brillouin zone with a technique developed in the 
lattice gauge theory^. 

The result for the Chern number in the present model 
obtained with the above method is displayed in Figl^l 
which is the key result in the present work, and re- 
veals an intriguing feature. For comparison, let us first 
look at Fig|2] (a), which is just the result for the ordi- 
nary graphene with the usual graphene Landau levels at 
= huj c sgii(n)y/n. The numerical result for the Chern 
number is simply the graphene QHE, where the well- 
known Chern number sequence of • • • , —3, —1, 1, 3, • • • 
for the honeycomb lattice coincides with a sum over the 
two, half-integer contributions from the degenerate pair 
of Dirac fields. 
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FIG. 2: (Color online) (a) For the ordinary graphene (ii = 
0) and (b) for the present model with a shifted Dirac cones 
(ti 7^ 0) we display the numerically calculated Chern numbers 
against Ef (upper panels), along with the Landau levels for 
each of the two Dirac cones, where the Chern numbers are 
indicated for Ef in each of the gaps (lower panels). The 
result is for the flux <j> = 1/100. 



If we now turn to the result for the present lat- 
tice model in Figj2] (b), the crucial question is: does 
the result conform to the expectation from the effective 
Dirac model where the sum of two half-odd-integer series 
(• • • , -3/2, -1/2, 1/2, 3/2, • • • ) with a shift in energy be- 
tween them should give the total Hall conductivity as 
E F traverses shifted sets of Landau levels (as displayed 
in the lower panel in the figure) . The numerical result for 
the Chern number in the lattice model is seen to agree 
with this. Namely, in a striking contrast to the ordinary 
graphene where each QHE step has a jump of 2 in the 
Chern number, the present model exhibits a jump of 1 at 
each step. This, along with the positions of these jumps, 
perfectly fit with the positions and the associated Chern 
numbers of the shifted Landau levels as the figure indi- 
cates. The agreement is rather surprising, since there is 
no obvious reason why the superposition of effective field 
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FIG. 3: (Color online) Numerical result (black curves) for the 
Landau level spectrum when the shift (A oc t\\ the horizon- 
tal axis) is varied, with the Chern number indicated for each 
of the Landau gap. The shading represents the Chern num- 
bers expected from the effective theory (see text), where the 
boundaries between different shades (Chern numbers) almost 
exactly coincide with the numerical result. Accumulation of 
lines around E = ±1 corresponds to the Van Hove singular- 
ity, which shifts quadratically with t\. The result is for a 
magnetic field of (f> = 1/100, and two vertical lines indicate, 
respectively, the situation for Figs[2ja) and[2jb). 



theory for the vicinities of K and K' and the lattice model 
should have the same Chern numbers. Thus, although we 
have still no half integers for the total Hall conductivity 
(since a sum of two half-odd integers is an integer), we 
have indirectly confirmed the half integers. Since no half 
integers can possibly appear according to TKNN, this is 
indeed as best as we can confirm the half-integer property 
of each Dirac cone. 

We expect such a decomposition of the Hall conduc- 
tivity into the contributions from each Dirac cone to sys- 
tematically occur when we vary the energy separation 
(2A) between the two Dirac points. FigJ3] displays the 
loci of Landau levels (solid lines) and the Chern number 
for each gap between the Landau levels when t\ oc A 
is varied. The result is just what is expected from the 
above decomposition, as indicated by the result for the 
two Dirac fields (shaded regions), which almost exactly 
coincides with the result for the lattice model. 
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IV. EDGE STATES 

The QBE state is a topological state, which imme- 
diately implies that edge states should appear at the 
boundaries of a sample, where an intimate relation ex- 
ists between the edge topological number (i.e., the edge 
Hall conductivity) and the bulk topological number (the 
bulk Hall conductivity)^. So it is an intriguing as well 
as necessary test to look into how this property appears 
in the edge-state in the present model. FigfJ] displays a 
numerical result for the energy spectrum against ki for a 
finite sample with zig-zag edges along X\. The spectrum 
comprises the bulk Landau levels (horizontal lines) and 
a series of edge modes (identified as curves that traverse 
between Landau levels). 

According to the bulk-edge correspondence in a topo- 
logical argument^, the bulk Chern number should coin- 
cide with the number of pairs of (right and left) edge 
modes in the Landau gap considered. In Fig|4] we can 
clearly recognize that the bulk Dirac Landau levels (oc 
y/n with the shift) (horizontal lines) are accompanied 
by right and left edge modes (upturning for the electron 
branch, down-turning for the hole branch). This occurs 
around each of K and K'. If we count the number of 
pairs of edge modes, we have, for each of the series, the 
beautiful sequence 1, 2, 3, • • • (as opposed to the sequence 
1, 3, 5, • • ■ for graphene, see Fig. 10 of—) both in the elec- 
tron or hole branches. The two series of pairs of edge 
modes arise with the energy shift between K and K' as 
they should. We can thus reinforce our argument, via 
the bulk-edge correspondence, that we can decompose 
the contribution from each of the Dirac cone that carries 
the bulk Chern number of 1/2 (Fig©. 

V. SUMMARY 

We have demonstrated that each Dirac 
fermion contributes the half odd integer series 
(••• ,-3/2,-1/2, 1/2,3/2,- ••) to the Hall conduc- 
tivity when the Fermi energy traverses the Landau levels 
in a lattice model in which the energy between the two 
Dirac points is systematically shifted. The picture is 
also confirmed, via the bulk-edge correspondence, in the 
edge mode of the Dirac fermions. 

Recently it is recognized that massless Dirac fermions 
can be realized as surface states of the three-dimensional 
topological (quantum spin Hall) systems^. There, the 
two-dimensional surface of the gapped three-dimensional 



topological insulator accommodates gapless (massless 
Dirac) fermions, which is a manifestation of the bulk-edge 
correspondence^. The massless Dirac fermions are again 
doubled for the topological consistency, but the doubling 
partner exists at the other side of the system in this case. 
Thus, whether we can decompose the topological num- 
bers into contributions from each Dirac cone is becoming 
a realistic question. It is an interesting future problem 
to consider how the present toy model can be realized 




FIG. 4: (Color online) Energy spectrum against fc2 in 
the present model for a finite sample with zigzag-bearded 
edge a 17 ' 18 . Horizontal solid (dashed) lines correspond to bulk 
Landau levels corresponding to the shifted Dirac fermions 
by A (—A). Curves are the edge modes whose amplitudes 
are confirmed to be localized along the zigzag right (red) or 
bearded left (blue) edge. The result is for a magnetic field 
<j> = 1/48 with ii/to = 0.02. 



in actual materials, where the surfaces of 3D topological 
insulators may be one possible avenue. Another possible 
realization of our model is a cleverly arranged cold atoms 
for an optical lattice where non-Abelian gauge structure 
can arise from the intra-atomic degree of freedom with a 
fine-tuning of desired parameters^. 
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